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ПАРАМЕТРАМИ В ЗАДАЧАХ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ 
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Рассматривается моделирование управления нестационарными температурными полями объектов со сложной геометрической 
конфигурацией. Приводятся сравнительные динамические характеристики замкнутых систем автоматического регулирования с 
перестраиваемой и фиксированной структурами. 


Введение 

Технические системы управления обычно явля- 
ются довольно сложными устройствами, динамика 
которых описывается различными функциональ- 
ными уравнениями. В каждом конкретном случае 
при использовании тех или иных математических 
методов необходимо составить математическую 
модель объекта. На практике подавляющее боль- 
шинство объектов - это объекты с распределенны- 
ми параметрами. Управляемый процесс с распре- 
деленными параметрами описывается краевыми 
задачами для дифференциальных или интегро- 
дифференциальных уравнений с частными произ- 
водными или бесконечными системами обыкно- 
венных дифференциальных уравнений [ 1] . 

Если рассматривать пространственно-много- 
мерные объекты со сложной формой границы обла- 
сти изменения пространственных координат, а так- 
же учитывать принципиально нелинейные эффек- 
ты, получить аналитическое решение уравнения 
объекта затруднительно. Данный аспект привел к 
широкому распространению на практике прибли- 
женных моделей объектов с распределенными пара- 
метрами упрощенного вида, описывающих их пове- 
дение с требуемой точностью. В инженерной прак- 
тике получили широкое распространение разност- 
ные методы приближенного описания объектов с 
распределенными параметрами, использующие раз- 
личные способы пространственного, временного 
или пространственно-временного квантования в 
области изменения аргументов входа и выхода рас- 
сматриваемого распределенного блока [2, 3]. 


Для построения моделей будем использовать ме- 
тод конечных элементов, позволяющий произвести 
дискретизацию области изменения пространствен- 
ных переменных путем разбиения с некоторой по- 
грешностью на ряд неперекрывающихся подобластей 
простой формы, в пределах каждого из которых функ- 
ция состояния объекта приближенно описывается од- 
нотипным образом линейной комбинацией конечно- 
го числа заранее выбранных базисных функций. 

Поддержание необходимого физического пара- 
метра на заданном уровне в подобных объектах яв- 
ляется важной и трудоемкой задачей. Одним из 
возможных способов обеспечения данного режима 
является вычислительная среда с перестраиваемой 
структурой, формирующая управляющее воздей- 
ствие, построенная по автоматному принципу и 
имеющая в узловых точках многофункциональные 
логические модули (МЛМ), т. е. автоматы с пере- 
страиваемой структурой [4]. 

Общее выражение для двувходового автомата 
можно представить в виде 

V = (АхА\<2,В,(р,у/). 

Функционирование такого автомата обуславли- 
вается некоторым сверхсловом ос ', ос ^(Л) 00 , и пред- 
полагается, что на первый вход автомата V (алфавит 
А) в каждый момент времени поступает произволь- 
ный входной сигнал, а на второй вход - только оче- 
редной символ сверхслова ос'. Следовательно, 
сверхслово ос 'управляет изменением структуры ав- 
томата V Под функционированием автомата с пе- 
рестраиваемой структурой V понимается тернарное 
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отношение {(а\(р(д,аха),у/(д,аха))\де (),аеА*}, 
где аха' - слово вида (а( 1),я\1))...(я(д),я'(д)), 
а=а(1)...а(8),а=а{1)а{2)... 

Поэтому разумное сочетание инженерной ин- 
туиции и математической строгости в анализе яв- 
ления составляет основу обоснованного выбора 
математической модели. 

Постановка задачи 

Исследуемый объект (рис. 1) состоит из трех ча- 
стей, каждая из которых занимает некоторую 
область Д (/-1,2,3 - номер части) трехмерного 
пространства, а Д - поверхность, ограничивающая 
соответствующую область Д, 


гм л 



Рис. 1. Исследуемый объект 


щее воздействие с ограничением: 0<//(/)<40; Д - 
поверхность, получаемая из выражения 


ь,=(с,\е')\еД 


где <21=(6,ПС І+І ). 




Рис. 2. Общий вид составных частей объекта: а) первая; 
б) третья ; в) вторая 


Функция Ні(х, у, г, і) описывает распределение 


температуры в соответствующей /-и части и удо- 
влетворяет уравнению теплопроводности 



где а { - коэффициент температуропроводности; х, 
у, г - пространственные декартовые координаты; 
р ; - плотность; с р . - теплоемкость; Я, - коэффици- 
ент теплопроводности. 

Для уравнения (1) начальное условие 

0 )= 0 , ( 2 ) 


граничные условия 



( 3 ) 

( 4 ) 


к 


і 


в! 




Здесь дк/дпУ означает производную по норма- 
ли п к поверхности 0/, взятую в точке, лежащей на 
поверхности ()} (рис. 2); / — номер части объекта; 
- коэффициент теплоотдачи; и(і) - управляю- 


Требуется, чтобы 

| к* -к(х\у\г\ оо) |< ^ оп , (6) 

где к* - требуемое значение температуры; х\у\і - 
заданные координаты точки А в трехмерном про- 
странстве (рис. 1); ДГ - допустимое значение ста- 
тической ошибки регулирования. 

Классический алгоритм управления 

Существует большое число способов [5] нахож- 
дения решения уравнения теплопроводности. В 
нашем случае мы имеем дело с уравнением, описы- 
вающим поведение геометрически сложного 
объекта в трехмерной декартовой системе коорди- 
нат. Так как для данного рода систем регулирую- 
щим органом является параметр объекта, в данном 
случае изменяющаяся температура на границах, то 
теоретическое исследование работы такого рода 
систем автоматического регулирования, особенно 
при наличии распределенных параметров, натал- 
кивается на большие математические трудности, 
связанные с исследованием нелинейных уравне- 
ний. Пожалуй, одним из эффективных способов 
достижения поставленной цели в сложившейся си- 
туации является привлечение численных методов и 
проведение имитационного моделирования. 
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Анализ и синтез современных систем автомати- 
ческого управления в настоящее время немыслимы 
без применения средств вычислительной техники. 
Существует большое количество пакетов приклад- 
ных программ (ППП), в которых реализованы по- 
пулярные и эффективные численные методы. Эти 
ППП способны облегчить решение задач синтеза и 
анализа подобных систем. К наиболее известным 
зарубежным ППП можно отнести МаіІаЬ, Майтсасі, 
Сотвоі МиШрЬуысз (РетЬаЬ), Маріе, Ріиепі. 

Наиболее интересным, по нашему мнению, яв- 
ляется пакет Сош8о1 МиШр1іу8ІС8, в основу которо- 
го заложен метод конечных элементов. Используя 
возможности ППП, произведем дискретизацию 
исходной геометрической модели (рис. 1), разби- 
вая ее на ряд конечных элементов. 

Сформируем управляющее воздействие и{і) по 
пропорциональному закону с ограничением: 


тах 


при ке{і)>к 


тах 5 


и(і) = \ к в (I) при к ■ <кв(і) < к 


' тах 5 


к тш при к е (і)< к - , 


где к шх , к т1п - максимально и минимально возмож- 
ное значение на выходе управляющего устройства 
соответственно; к — коэффициент передачи упра- 
вляющего устройства; в(і) - ошибка регулирования. 




Рис. 3. Температурное поле на момент времени: а) 300 и 
6) 655 с 


Зададим к= 16; к ш = 40; к тт =0; а { =ау= 7,69- ІО 6 ; 
я 2 =8,9*10“ 7 ; Н*= 10; г с д т оп =0,5; Я 1 =Я 3 =27; Я 2 =1,3. Здесь 
и далее, когда единицы измерения не указаны, то 
числа представляют соответствующие физические 
величины в любой согласованной системе единиц. 

В результате проведения имитационного моде- 
лирования получено распределение температурно- 
го поля по поверхности исследуемого тела (рис. 3). 

Динамика изменения температуры в контроль- 
ной точке (х\у\і) трехмерного пространства пред- 
ставлена на рис. 4, а , а управляющее воздействие, 
при котором получена эта динамика, - на рис. 4, б. 

Из рис. 4, а видно, что в системе присутствует 
статическая ошибка ^ ст =0,5, возникающая в ре- 
зультате использования пропорционального зако- 
на регулирования в статической системе и удовле- 
творяющая заданному требованию (6). Чтобы све- 
сти ее к нулю, достаточно в контур управления до- 
бавить интегрирующее звено, тем самым наделив 
всю систему свойством астатизма. 




Рис. 4. Динамические характеристики: а) переходный про- 
цесс в контрольной точке; б) управляющее воздей- 
ствие 

Итак, в системах с распределенными параметра- 
ми допустимо применять законы управления, 
функционирующие в системах с сосредоточенными 
параметрами, в случае, если предъявляемые требо- 
вания к системе ограничиваются поддержанием 
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технологического параметра в одной единственной 
точке и никаких других ограничений на соседние 
точки технологического объекта не наложено. 


Структурно-перестраиваемый алгоритм управления 


Рассмотрим случай, когда динамика распро- 
странения тепла в исследуемом объекте описывает- 
ся выражениями (1)-(4) с дополнительным гра- 
ничным условием 



= п 1 (н(1)~М л < ), 


где д(1) - выходной сигнал исполнительного меха- 
низма, описываемый уравнением 


Т им ц\і) = и(1). 


Здесь Т им - постоянная времени исполнитель- 
ного механизма; и(1) - управляющее воздействие. 

Обеспечим в системе переходный процесс с за- 
данным качеством (время регулирования 60 с, пе- 
ререгулирование 20 %) в точке с координатами х\ 
у \ I. Для этого сформируем управляющее воздей- 
ствие, используя в управляющем устройстве, бази- 
рующемся на МЛМ, принцип перестраиваемости 
структур [6, 7]: 

и(1 ) = Ч* 

Г а при в(1)8 > 0, 

\р при 8(1)8 < 0, 


где а и /3- коэффициенты передачи первой и вто- 
рой линейных структур, используемых в управляю- 
щем устройстве соответственно; 8=$%п(8'(1)+С8(1)) 
- информация о знаке линейной комбинации 
ошибки и ее производной, характеризующая поло- 
жение системы в фазовом пространстве относи- 
тельно прямой переключения; с - коэффициент 
наклона прямой. Как видно, для функционирова- 
ния системы не требуется точного значения произ- 
водной от сигнала ошибки, а достаточно лишь ин- 
формации о знаке ее линейной комбинации с ве- 
личиной ошибки, которую можно получить срав- 
нительно простыми техническими средствами, на- 
пример, описанными в [8]. 

Зададим значение коэффициента а максималь- 
ным, которое может обеспечить управляющее 
устройство, а коэффициент (3 достаточно выбрать 
так, чтобы система в момент переключения на вто- 
рую линейную структуру была охвачена положи- 
тельной обратной связью и ее фазовая траектория 
обращена в противоположное направление движе- 
ния траектории первой линейной структуры (при 
коэффициенте а). Итак, а=50; Д=-10; с=0,1; Т им = 5. 

Из результатов моделирования, рис. 5, следует, 
что полученный переходный процесс удовлетворя- 
ет заданным оценкам качества: время регулирова- 
ния составляет 58 с и перерегулирование - 0 %. 

При рассмотрении первой системы автоматиче- 
ского регулирования с классическим алгоритмом 
управлением (на технических средствах Іпіеі Репіі- 


шп М 1,7 ГГц, ОЗУ 2 Гб) моделирование динамиче- 
ских характеристик длительностью 900 с заняло 10 
мин, а для системы со структурно-перестраивае- 
мым алгоритмом управления 90 с - 3,5 ч. Связано 
это с тем, что в структурно-перестраиваемой систе- 
ме в граничные условия, описывающие динамиче- 
ские характеристики исследуемого объекта, входит 
разрывная функция, т. е. динамика процесса стано- 
вится принципиально нелинейной. Следователь- 
но, при имитационном моделировании подобных 
систем накладываются более жесткие требования 
на аппаратные средства, с помощью которых про- 
водится вычислительный эксперимент. 



Рис. 5. Динамические характеристики в системе с пере- 
страиваемой структурой: а) переходный процесс; 
6) регулирующее воздействие на объект 


Заключение 

Приведены результаты математического моде- 
лирования нестационарных тепловых полей регу- 
лируемых объектов сложной геометрии, когда в 
граничные условия входит функция регулирования. 

На основе полученной математической модели 
проведены вычислительные эксперименты по ис- 
следованию динамики систем автоматического ре- 
гулирования с двумя типами алгоритмов управле- 
ния: классическим (пропорциональным) и струк- 
турно-перестраиваемым. В результате исследова- 
ний выявлено, что моделирование последних зани- 
мает значительно больше времени вычислений, 
чем в классических системах. Однако проигрыш во 
времени моделирования при анализе системы ком- 
пенсируется значительным выигрышем в динами- 
ческих характеристиках функционирующей систе- 
мы (быстродействие в 10 раз выше и отсутствует 
перерегулирование) . 
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Применение декартовой системы координат в по- 
строенной модели позволяет производить исследова- 
ния объектов как с осью симметрии, так и без нее. 
Особое внимание при рассмотрении подобных задач 
следует уделять построению геометрии исследуемого 
объекта, от которой зависит как скорость вычисле- 
ния, так и вообще их возможность. Следует избегать 
избыточности мелких деталей, оказывающих малое 
влияние на динамику изучаемого процесса. 
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